Observing the dark matter density profile of isolated galaxies 
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ABSTRACT 

Using the Sloan Digital Sky Survey (SDSS), we probe the halo mass distribution by studying 
the velocities of satellites orbiting isolated galaxies. In a subsample that covers 2500 sq. degrees on 
the sky, we detect about 3000 satellites with absolute blue magnitudes going down to Mb = — 14; 
most of the satellites have Mb = — 16 to —18, comparable to the magnitudes of M32 and the 
Magellanic Clouds. After a careful, model-independent removal of interlopers, we find that the 
line-of-sight velocity dispersion of satellites declines with distance to the primary. For an 
galaxy the r.m.s. line-of-sight velocity changes from «120 km s _1 at 20 kpc to «60 km s _1 
at 350 kpc. This decline agrees remarkably well with theoretical expectations, as all modern 
cosmological models predict that the density of dark matter in the peripheral parts of galaxies 
declines as Pdm oc r -3 . Thus, for the first time we find direct observational evidence of the density 
decline predicted by cosmological models; we also note that this result contradicts alternative 
theories of gravity such as MOND. We also find that the velocity dispersion of satellites within 
100 kpc scales with the absolute magnitude of the central galaxy as a oc L 3 ; this is very close 
to the Tully-Fisher relation for normal spiral galaxies. 

Subject headings: galaxies:dwarfs — galaxies:halos — galaxies:kinematics and dynamics — dark matter 
— surveys 



1. Introduction 

Measuring the distribution of mass around 
galaxies at large radii (> 100 kpc) provides a 
critical test for cosmological models. By study- 
ing the mass distribution at these distances, we 
directly address one of the most intriguing ques- 



tions of modern cosmology - the nature of dark 
matter. Specifically, we wish to determine the ex- 
tent of galaxy dark matter halos and their density 
profile. While the inner parts of galaxy halos have 
density distributions that yield approximately flat 
rotation curves, our understanding of the profile 
at larger distances is much poorer. The issue 
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is critical because the density profile that gives 
rise to a flat rotation curve (p oc r~ 2 ) is differ- 
ent from that predicted by cosmological models 
(NFW Navarro, Frenk, & White 1997) at larger 
distances (p oc r~ 3 ). 

The challenge in measuring the mass distribu- 
tion at large radii arises from the difficulty in find- 
ing a visible tracer to probe the mass. Historically, 
neutral hydrogen (HI) has been used to study the 
outer parts of rotation curves (e.g., Roberts & 
Rots 1973; Bosma 1981; van Albada et al. 1985), 
since HI is detected well beyond the optical bound- 
aries of spiral galaxies. The observations demon- 
strating that HI rotation velocities of field spi- 
rals are roughly constant at large galactocentric 
distances provides one of the primary pieces of 
evidence that galaxies are embedded in massive 
dark matter halos (e.g., Ostriker, Peebles, & Yahil 
1974; Faber & Gallagher 1979; Sancisi & van Al- 
bada 1987; Sofue & Rubin 2001). However, HI 
emission is detected only out to 30-50 kpc. X-ray 
emission from the diffuse hot gas in isolated ellip- 
tical galaxies also indicates the presence of mas- 
sive dark matter halos (e.g., Buote & Canizares 
1994; Romanowsky & Kochanek 1998; Buote et al. 
2002), but, just like HI in spiral galaxies, the X- 
ray emission does not extend far enough to distin- 
guish isothermal from cosmological p oc r~ 3 pro- 
files (Buote et al. 2002). Strong gravitational lens- 
ing places an important upper limit on the amount 
of dark matter but only in the inner (few tens of 
kpc) regions of galaxies (e.g., Keeton, Kochanek, 
& Falco 1998; Keeton 2001; Mailer et al. 2000). 

Weak gravitational lensing provides a more 
promising method of studying the outer parts of 
galaxies. Individual field galaxies produce a small 
distortion of background galaxies, and this dis- 
tortion can be used to measure the galaxy-mass 
correlation function. Unfortunately, the existing 
data (e.g., Tyson et al. 1984; Brainerd, Bland- 
ford, & Smail 1996; Hoeskstra, Yee, & Gladders 
2001; Fischer et al. 2000; Smith et al. 2001) do not 
distinguish between a singular isothermal sphere 
(SIS; Schneider & Rix 1997) and a cosmological 
NFW profile: both models provide good fits to the 
galaxy-mass correlation function in the outer re- 
gions where the error bars are large. Moreover, the 
galaxy-mass correlation function can be affected 
by neighbors of the lens galaxy (e.g., Smith et 
al. 2001; Guzik & Seljak 2002). Once the con- 



tamination from the neighboring galaxies is taken 
into account, as done by, for example, Fischer et 
al. (2000); Hoeskstra, Yee, & Gladders (2001), 
McKay et al. (2001) place only a lower limit on 
the size of the halo. 

Velocities of satellites of galaxies provide an- 
other way to probe the mass distribution at large 
radii. These have been used to constrain the 
mass of the Milky Way (Zaritsky et al. 1989; Lin, 
Jones, & Klemola 1995; Kochanek 1996; Evans 
& Wilkinson 2000) and M31 (Evans et al. 2000). 
However, because the number of detectable satel- 
lites around galaxies outside of the Local Group 
is small, one needs to study many galaxies to ac- 
cumulate enough statistics to study the profiles of 
other galaxies. As a result, observational efforts to 
study the dynamics of satellites have been some- 
what limited (Erickson, Gottesman, & Hunter 
1987; Zaritsky et al. 1993; Zaritsky & White 1994; 
Zaritsky et al. 1997; McKay et al. 2002). Impor- 
tant early results were obtained by Zaritsky et al. 
(1993), Zaritsky & White (1994) and Zaritsky et 
al. (1997), who compiled and studied a sample of 
about 100 satellites of nearby isolated spiral galax- 
ies with an average of 1-2 satellites per primary 
galaxy. It was found that the line-of-sight veloc- 
ity dispersion of the satellites does not decline with 
the projected distance to the primary galaxy. This 
result has been generally considered a strong ar- 
gument for the presence of dark matter at large 
distances (~200-400 kpc) around galaxies. Zarit- 
sky & White (1994) also found that the satellite 
velocity dispersion does not correlate with the lu- 
minosity of the primary galaxy. 

McKay et al. (2002) used the SDSS to select 
and study a much larger sample of 1225 satel- 
lites. The mean number of satellites they found 
around each host was only two. The new analysis 
confirmed that the velocity dispersion of satellites 
does not decline with distance from the primary, 
as observed by Zaritsky et al. (1993); Zaritsky & 
White (1994); Zaritsky et al. (1997). This implies 
that halos of isolated galaxies extend to distances 
of several hundred kpc. In contradiction with 
the earlier results, however, McKay et al. (2002) 
find that the velocity dispersion of satellites, cr, 
increases with the luminosity, L, of the primary 
galaxy as cr oc L 5 and claim that this is con- 
sistent with the predictions of the semi-analytical 
models of galaxy formation of Kauffmann et al. 
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(1999) when the sample of primary galaxies in the 
models is defined as it is in the observational anal- 
ysis. 

Cosmological models made definite predictions 
for the mass profile of dark matter halos at large 
distances. Although SIS halo models with den- 
sity p cx r~ 2 and a = constant are often assumed, 
there is little justification for using SIS at these 
distances. The main motivations for using SIS are 
simplicity and extrapolation of flat rotation curves 
to much larger radii. However, during the last two 
decades of intensive numerical modeling of galaxy 
formation, not a single model has produced an 
SIS. Every cosmological model studied so far (and 
there have been plenty) has p cx r~ 3 at large radii. 
The slope does not depend on the mean density 
of matter: CDM with f2 ma tter = 1 (Ghigna et al. 
2000) and A CDM with f7 mattCT = 0.3 (Klypin et 
al. 2001) have the same slope. It does not depend 
on the nature of the dark matter: warm dark mat- 
ter (Avila- Reese et al. 2001), self-interacting dark 
matter (Colin et al. 2000), and cold dark matter 
all make the same prediction. The slope does not 
depend on the halo mass: halos ranging from clus- 
ter masses to dwarf masses all have p cx r~ 3 . In 
other words, the declining velocity dispersion of 
satellites is not a test for the parameters of cos- 
mological models, but rather of the hierarchical 
scenario itself. It is ironic that the same argument 
- the constant velocity dispersion of satellites - 
which just a few years ago provided one of the 
strongest cases for the existence of dark matter, is 
now an argument against dark matter. The only 
model that predicts a constant velocity dispersion 
is Modified Newtonian Dynamics (MOND). 

As previous authors have noted (Zaritsky 1992; 
Zaritsky & White 1994), a comprehensive mea- 
surement of the satellite velocity dispersion pro- 
file around isolated galaxies should ideally com- 
bine two main elements: (1) a large number of 
primary galaxy and satellite (candidates) to esti- 
mate the satellite kinematics in separate bins of 
primary mass (or luminosity) and projected sepa- 
ration, and (2) a model- independent, robust, and 
testable rejection of interloper galaxies - dwarfs 
with large physical distances from the primary 
galaxies but small projected and velocity differ- 
ences - that are not bound to the primary. The ef- 
fect of the interlopers is to make the halo mass pro- 
file difficult to measure, and that may have led to 



systematic effects in the few observations that sug- 
gest a constant velocity dispersion at large radii. 
Understanding their effect is crucial for any inter- 
pretation of the halo mass distribution. Results 
from weak lensing can be affected by the presence 
of nearby projected neighbors, leading to similar 
systematic problems (see Smith et al. 2001; Guzik 
& Seljak 2002). 

In this paper we use the vast data base of SDSS 
to study the motion of satellites around a care- 
fully selected sample of isolated nearby galaxies. 
We also develop an improved, model-independent 
approach for interloper rejection, which, together 
with state-of-the art cosmological iV-body simula- 
tions of galaxies and their satellites, allows us to 
measure the dark matter halo profile. 

In Section 2, we describe the observational data 
and the criteria for selecting our samples of pri- 
maries and satellites. The numerical simulations 
are briefly described in Section 3, where we also 
discuss the use of simulations to test prescriptions 
to remove interlopers. Our observational results 
are presented in Section 4, and our conclusions 
are given in Section 5. 

2. Observational Data: Selection of pri- 
maries and satellites 

The SDSS (York et al. 2000; Stoughton et al. 
2002) is a survey that images up to 10 4 deg 2 of 
the the northern Galactic cap in five bands ugriz 
(Fukugita et. al. 1996; Pier et al. 2002; Hogg et. 
al. 2001; Smith et al. 2002) using a drift-scanning 
mosaic CCD camera (Gunn et al. 1998), down to a 
limiting magnitude of r' <~ 23 (Lupton et al. 2001). 
Approximately 900,000 galaxies down to r' ~ 17.7 
will be targeted for spectroscopic follow-up using 
two fiber-fed spectrographs on the same telescope 
(Strauss et al. 2002; Blanton et al. 2002). 

For the current analysis, we use data on galax- 
ies for which spectra were obtained before Au- 
gust 2002. We use "survey quality" redshifts de- 
rived from these spectra by D. Schlcgcl. The 
recessional velocity errors are always less than 
20 km s _1 . Our SDSS rcdshift sample consists 
of 254,073 galaxies distributed in several strips 
on the sky; the total area covered by our data is 
about 2,500 deg 2 . Photometric parameters were 
taken from the PHOTO measurement (Lupton et 
al. 2001) of the brightest object within 3 arc- 
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Fig. 1. — Sky coverage for the SDSS dataset used 
in this work. The total area covered by our data 
is -2500 deg 2 . 

sec that was loaded into the SDSS collaboration 
database in August 2002 (but see also discus- 
sion below); magnitudes were corrected for fore- 
ground extinction using the (Schlegel, Finkbeiner, 
& Davis 1998) values loaded in the SDSS database. 
The sky coverage of the resulting SDSS imaging 
and spectroscopic data is shown in Figure 1. 

We supplemented the SDSS data with data 
from the RC3 catalog (de Vaucouleurs et al. 1991) 
for two reasons. First, as we are identifying iso- 
lated galaxies (see below) , we want to avoid choos- 
ing galaxies that have bright neighbors outside the 
area covered by the SDSS at the time we made 
our sample; the RC3 provides a catalog of galaxies 
over the whole sky to make this possible. While 
it is true that the RC3 catalog is not complete, 
especially at larger distances, missing companions 
in the SDSS data are much more likely for the 
nearest galaxies for which the area on the sky cor- 
responding to our isolation criterion is large. For 
these, the RC3 catalog should be nearly complete. 
Another, possibly more important, reason for us- 
ing the RC3 is that the version of PHOTO that 
we were using did not provide accurate photom- 
etry for large and bright galaxies owing to prob- 
lems with the deblending algorithm. As a result, 
we used RC3 Bt magnitudes (or photographical 
magnitudes if Bt was not available) , corrected for 
foreground reddening, for all of the SDSS galaxies 
for which we found a match in the RC3 catalog 
(within 0.5 arcmin on the sky and 100 km s^ 1 in 
velocity); for all of the other galaxies (which are 
the smaller galaxies for which the SDSS photome- 



try is accurate), we transformed the SDSS g mag- 
nitude to B magnitudes using the relation given 
by Yasuda et al. (2001). 

Our main galaxy list was selected from the 
full redshift sample by taking all galaxies with 
500 < cz < 60,000 km s _1 and apparent magni- 
tude r' < 17.7; the total number of selected galax- 
ies is 206,352. Distances are derived assuming a 
smooth Hubble flow with h = 0.7 (H = 100 hr Y 
km s _1 Mpc -1 ); the heliocentric velocities were 
converted to the Local Group Standard of Rest 
before computing distances. 

For the study of satellite dynamics, we de- 
fine primary galaxies to be those with absolute 
blue magnitude brighter than Mb — —19.0. We 
put no restriction on Hubble type. We use three 
different isolation and selection criteria to define 
three different samples of isolated primaries. To 
be considered as isolated, a galaxy must have no 
other galaxies within a magnitude difference AM, 
projected separation AR, and velocity separation 
AV. Satellites are defined as all objects within 
a projected distance Sr and velocity difference 6v 
and being at least 5m magnitudes fainter than the 
primary. Table 1 gives the parameters used to 
define each sample, and the number and median 
distance of primaries, and the number of satellites. 

The first two samples have identical isolation 
criteria, but Sample 2 is defined for primaries out 
to a larger distance, hence only bright satellites 
can be seen for the more distant galaxies. The 
shallower sample has more satellites per primary, 
but the more distant sample has an overall larger 
number of total satellites. The third sample mim- 
ics conditions used by McKay et al. (2002), in 
which the isolation criterion is significantly re- 
laxed. It has a very large search radius (almost 
3 Mpc), but the luminosity of the neighbors can 
be quite large: half that of the luminosity of the 
primary. In other words, this prescription can pick 
a bright member of a group of galaxies and erro- 
neously treat it as an "isolated" galaxy. As dis- 
cussed below we obtain similar results for all of 
the samples, suggesting that our results are statis- 
tically significant and do not depend on the exact 
isolation criterion. Our main results are based on 
the first two samples. 

Figure 2 shows distributions of primary and 
satellite brightnesses, brightness ratios, and num- 
ber of satellites per primary for Sample 1. Consid- 
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Table 1 

Selection and isolation criteria for SDSS samples 



Parameter 


Sample 1 


Sample 2 


Sample 3 


Alaximum depth of the sample (km s ) 


1U,UUU 


fin nnn 


fin nnn 
ou,uuu 


Constraints on bright neighbors; 








Magnitude difference, AM 


z.U 


o n 

z.U 


n 7^ 

U. I 


Minimum projected distance, Ai? (/i 1 kpc) 
Minimum velocity separation, Ay (km s _1 ) 


^nn 
ouu 


ouu 


onnn 
zUUU 


1000 


1000 


1000 


Constraints on satellites: 








Minimum magnitude difference, dm 


2.0 


2.0 


1.5 


Maximum projected distance to primary, Sr (h^ 1 kpc) 


350 


350 


500 


Maximum velocity separation with primary, Sv (km s _1 ) 


500 


500 


1000 


Number of isolated galaxies 


1278 


88603 


26807 


Number of satellites 


453 


1052 


2734 


Statistics of isolated galaxies with at least one satellite: 








Number 


283 


716 


1107 


Mean distance (km s _1 ) 


7100 


14700 


23170 


Mean distance (km s" 1 ) for -19.5 < M B < -20.5 
Mean distance (km s" 1 ) for -20.5 < M B < -21.5 


7244 


9785 


11076 


7697 


15917 


20854 


Limiting magnitude Mb 


-16.4 


-18.0 


-19.0 



cring only primaries for which there is at least one 
satellite, we find on average two satellites per pri- 
mary in our first sample, which is similar to that 
found by Zaritsky et al. (1997) and McKay et al. 
(2002). A visual morphology classification of our 
primaries in Sample 1 indicates that about 30% of 
the galaxies are ellipticals and SOs (E-SO). This 
fraction of early-type galaxies is simiar as seen in 
the RC3 catalog and in the SDSS (Nakamura et 
al. 2002). A fraction of our primaries is close to 
the edges of the strips and therefore we do not 
cover their entire halos; for 84% of the primaries 
the whole volume is covered. 

Given the limiting magnitude of the SDSS spec- 
troscopic survey, we expect to find satellites down 
to Mb ~ —14 for a primary with cz = 2000 
km s _1 and down to Mb ~ —17 for a primary 
with cz = 7100 km s _1 , the mean of our pri- 
mary recessional velocities in Sample 1. These 
limits would yield about four satellites per pri- 
mary for our nearest primaries and 1-2 satellites 
for the further ones if we assume that all galaxies 
have satellite luminosity functions similar to that 
of the Local Group (see Pritchet, & van den Bergh 



1999; Grebel 2000). Samples 2 and 3 include more 
distant primaries, and only brighter satellites for 
the more distant galaxies (see Table 1 for limiting 
magnitudes) . 

The main physical properties of satellites and 
primaries, such as their color distribution, spectral 
properties, and their subdivision into early- and 
late-type galaxies will be discussed in Vitvitska et 
al. (2003). 

3. Numerical models 

Cosmological simulations of galactic dark mat- 
ter satellites in the hierarchical model of struc- 
ture formation predict a remarkably large number 
of dark matter satellites orbiting around a Milky 
Way-size halo (Klypin et al. 1999; Moore et al. 
1999). In many respects the simulations are quite 
realistic; we use them to develop a method for re- 
moving the effects of interlopers on the amplitude 
and shape of the r.m.s. velocity profile of the satel- 
lites. Note that the results of the simulations are 
not used in the reduction of the real astronomical 
data. 

We use one of the simulations presented by 
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AM B (satellite-primary) Number of satellites per primary 

Fig. 2. — Statistics of primary galaxies and their 
satellites in a subsample with velocities less than 
10, 000 km s _1 . The top left panel shows the dis- 
tribution of absolute blue magnitudes, Mb, of iso- 
lated primaries that have at least one satellite; the 
top right panel shows the brightnesses of the satel- 
lites. The magnitude differences between the pri- 
maries and the satellites is presented in the bottom 
left panel. The number of satellites around each 
host is shown in the bottom right panel. 

Klypin et al. (2001). The simulation is done for 
the standard A CDM cosmological model with 
ilo = 0.3, Cl\ = 0.7, and h = 0.7; the spec- 
trum normalization is cr 8 = 0.9. The Adaptive- 
Rcfinement-Tree (ART) code (see Klypin et al. 
2001, for details) is used to run the simulation. 
The simulation box was 25 h~ x Mpc. The formal 
force resolution was 100 /i -1 pc; the mass resolu- 
tion is 1.2 x 10 6 h- 1 M Q . 

The simulation has three galaxy-size halos. All 
halos are relaxed and have cusps. The outer parts 
of the density profiles of the halos are well approx- 
imated by an NFW profile. The distance between 
two of the halos is 600 kpc; we decided not to use 
this pair because the halos are not isolated. We 
selected the third halo, which does not have a mas- 
sive companion within 3 Mpc radius. This halo 
has a virial mass of 1.5 x 10 12 h^ 1 M Q , a virial 
radius of 235 h" 1 kpc, and a maximum circular 



velocity of 214 km s _1 . There are 200 dark mat- 
ter satellites inside a 350 h~ x kpc radius from the 
center of the halo. Most of the satellites are very 
small with circular velocity of 10-15 km s _1 and 
mass 5 x 10 7 h^ 1 M . 

The simulations were done in such a way that 
only regions around the three large halos ("~l/2 
Mpc) have high resolution. The rest of the compu- 
tational volume was simulated with significantly 
lower resolution. No halos outside the regions of 
high resolution are used in our analysis. 

3.1. Satellites in cosmological simulations: 
removal of interlopers 

In any sample of isolated bright galaxies and 
their satellites there is a fraction of galaxies that 
are probably not related to the host, but which are 
included in the sample because of projection ef- 
fects. We call these objects interlopers. The inter- 
lopers tend to have large velocity differences and 
large projected separations. The fraction of inter- 
lopers in the samples of Zaritsky & White (1994) 
and Zaritsky et al. (1997) was 15%. We expect 
a similar level of contamination in our observa- 
tional samples. If the interlopers are not removed, 
the r.m.s. velocity of satellites, and the inferred 
mass of the primaries, will be overestimated. We 
use the simulations described in the last section to 
test prescriptions for removing the interlopers. 

To understand how interlopers might affect the 
results we introduce them into the simulation. We 
have to do this because the A^-body simulation 
resolves only the region very close to the central 
halo; no interlopers are present in the simulation. 
To make a sample that realistically mimics the ob- 
servational data we start by making an estimate of 
the number of expected interlopers in the observed 
sample. Using the SDSS field luminosity function 
(Blanton et al. 2001) we find that 0.43 galaxies 
between — 18 < Mb < — 16 are expected inside 
a volume R < 350 hr x kpc and with a velocity 
dispersion of |A1^| < 500 km s _1 at a distance of 
^5000 km s _1 . This is the estimated number of 
interlopers per primary, if galaxies are randomly 
distributed in space. Neglecting the large-scale 
correlations of galaxies appears to be a reasonable 
assumption at this point, because our primaries 
are selected to avoid clusters and groups of galax- 
ies. The latter are responsible for most galaxy 
clustering on a scale of a few megaparsecs. Thus, 
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Fig. 3. — The absolute value of the line-of-sight 
velocity difference between satellites and the cen- 
tral halo as a function of projected distance to 
a simulated dark matter halo in the standard 
A CDM model. Filled circles show (true) satel- 
lites with 3D distance less than 350 hr 1 kpc from 
the central halo with mass M v ; r = 1.5 x 
10 12 ft- _1 Mq. Open circles show a realization 
of interlopers, which were generated by ran- 
domly placing objects inside a cylinder of radius 
350 hr 1 kpc and length 1000 km s" 1 . The num- 
ber of interlopers is defined by the luminosity func- 
tion of galaxies in the SDSS. The dashed curve 
gives the r.m.s. velocity in the combined sample. 
The full curve shows the r.m.s. velocity of true 
satellites. 

the clustering of galaxies (with the exception of 
true satellites) around our primaries is small and 
can be neglected. Using simulations this point 
should be quantified in more detail in the future. 

In the observed samples we find on average two 
satellites per primary. Thus, we expect 0.215 in- 
terlopers for each true satellite. In our numeri- 
cal simulation we have 200 true satellites. This 
means that we expect about 43 interlopers. We 
randomly place 43 interlopers into a cylinder of 
radius R = 350 hr 1 kpc and length 1000 km s _1 . 

Figure 3 shows the resulting velocity differences 
for both satellites and interlopers as a function of 



projected distance to the parent dark matter halo, 
as well as the r.m.s. velocities for the true satellites 
and for the combined sample of the satellites and 
interlopers. Since the interlopers are not physi- 
cally associated with the primaries, their number 
depends only on the volume. Therefore, their ef- 
fect on the r.m.s. velocity is stronger at larger 
distances. Neglecting interlopers clearly leads to 
an apparent increase in the r.m.s. velocity with 
distance, while the r.m.s. velocity of the true satel- 
lites actually decreases with distance. 

To reliably estimate the satellite velocity dis- 
persion, one needs to find a way to account for 
the interlopers. We start by characterizing the 
velocity distribution of true satellites. Motivated 
by Figure 3, we consider a Gaussian velocity dis- 
tribution, but allow the width of the distribution 
to change with radius. Figure 4 shows the re- 
sulting distribution of normalized (to the r.m.s. 
velocity in different radial bins) velocities, and 
demonstrates that the distribution is very close 
to Gaussian. In reality, the distribution cannot be 
completely Gaussian because satellites with very 
high velocities will not be bound to the halo, and 
will escape. Using simple arguments based on the 
virial theorem, one naively expects that satellites 
with velocities twice the r.m.s. velocity should not 
be bound. Indeed, we find an indication of this: in 
Figure 4 the points are systematically below the 
Gaussian fit for AV > 2a. However, it is difficult 
to use the simple argument because we measure 
only one component of the velocity and because 
the escape velocity changes with distance. In any 
case, the Gaussian fit provides a reasonably ac- 
curate approximation for the distribution of the 
line-of-sight velocities of true satellites. 

A simple prescription far removing the interlop- 
ers was used by McKay et al. (2002). The veloc- 
ity difference distribution of the entire sample of 
neighbors (including objects at all projected dis- 
tances) was fitted by a Gaussian + constant func- 
tion. The idea is that the Gaussian describes the 
velocity distribution of satellites, and that the con- 
stant takes into account the distribution of inter- 
lopers; the width of the Gaussian component gives 
a direct measurement of the r.m.s. velocity of the 
satellites. Our simulations show that the velocities 
of the true satellites are indeed well approximated 
by the Gaussian distribution, but that the width 
of the distribution (i.e., the velocity dispersion) is 
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Fig. 4. — The distribution of relative velocities of 
true satellites in TV-body simulations. The distri- 
bution is very close to Gaussian (full curve). The 
distribution was obtained by finding the r.m.s. ve- 
locities for individual bins and then combining the 
normalized velocities. 

a function of radius. We find that the McKay et al. 
(2002) prescription has two problems that result 
from making a single fit over all radii: 

(1) The level of interloper contamination (the 
constant in the fit) does not depend on the dis- 
tance to the primary. This is inaccurate because 
of a purely geometrical factor: the volume and 
thus, the number of interlopers, is larger for larger 
radii. Neglecting this dependence produces a bias: 
the r.m.s. velocity is overestimated at large radii 
and underestimated at small radii. 

(2) There is no allowance for dependence of the 
velocity dispersion on projected distance. Because 
this dependence is in fact what we wish to study, 
we do not want to impose the assumption of a 
constant velocity dispersion on the prescription for 
the removal of interlopers. 

In fact, if we apply the McKay et al. (2002) pre- 
scription to the simulations, we find that the re- 
covered profile of the (interloper-removed) r.m.s. 
velocity is flat; the prescription fails to recover the 
true dispersion profile of the satellites. In this ar- 
ticle, we apply an improved interloper prescrip- 



tion by making a separate "Gauss + constant" fit 
for separate radial bins. In this case the increase 
of the number of interlopers with the projected 
distance is taken into account (at least approxi- 
mately) . This also allows for a radially dependent 
velocity dispersion of true satellites. In practice we 
use 100 kpc radial bins (70 h^ 1 kpc in numerical 
simulations) and 100 km s _1 velocity bins. We do 
not include satellites with projected distance less 
than 20 kpc. The binned data are used to mini- 
mize the sum over all bins of normalized deviations 
(Vi,obs - ^modei) 2 /^, where iV; is the number of 
velocities in the zth bin, and l^, bs and Vi. mo dei are 
the observed and modeled velocities in the bin. 

Figure 5 shows the accuracy of the recovery 
of the true r.m.s. velocity for the simulation; 
the short-dashed line shows the recovered veloc- 
ity while the solid line shows the true velocity dis- 
tribution of the satellites. The procedure works 
remarkably well. It recovers the value of the ve- 
locity dispersion and the number of interlopers. In 
the next section, we apply the same procedure to 
the observational samples. 

The theoretical prediction for the velocity dis- 
persion profile for an NFW profile with the same 
virial mass as that of the simulated dark matter 
halo is shown by a long-dashed line. In order to 
make this prediction for the r.m.s. line-of-sight 
velocity, we use the NFW profile and assume that 
the satellites are in equilibrium. We then use the 
Jeans equation to find the velocity dispersion. The 
velocity anisotropy (3(r) = 1 — a]_/2<j^., which is 
needed for the equation, was taken from cosmolog- 
ical simulations (Colin, Klypin, & Kravtsov 2000; 
Vitvitska et al. 2002). (Here a± and ay are the 
tangential and radial r.m.s. velocities). The ve- 
locity anisotropy changes from very small values 
(3 w close to the center of a halo to j3 ~ 0.5 — 0.6 
at the virial radius. We use this dependence of 
(3 on radius to make analytical estimates of the 
projected velocity dispersion. Specifically, we in- 
tegrate the r.m.s. velocities, which we obtain from 
the Jeans equation, along a line of sight. The in- 
tegration is truncated at 1.5 virial radii. The ve- 
locity anisotropy enters these calculations twice: 
through the Jeans equation and during the inte- 
gration along the line of sight. For a fixed mass 
profile M(r) the effects of the velocity anisotropy 
are rather mild. The line-of-sight r.m.s. velocity 
in the realistic non-isotropic halo declines with dis- 
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Fig. 5. — The same as in Figure 3. The long 
dashed and dot-dashed curves show respectively 
the velocity dispersion and the escape velocity for 
the NFW profile with the same virial mass as the 
real dark matter halo. Most of the interlopers are 
above the escape velocity - they cannot possibly 
belong to the dark matter halo. The r.m.s. veloc- 
ity of the true satellites is plotted using the solid 
curve. This r.m.s. velocity of true satellites is ac- 
curately recovered (short-dashed curve) from the 
total sample of satellites + interlopers by making 
the Gaussian + constant fit for velocities in indi- 
vidual radial bins. 

tance slightly more rapidly as compared with the 
isotropic case (/? = 0). 

4. Results 

4.1. Radial dependence of satellite veloc- 
ity dispersion 

Figure 6 shows the relative line-of-sight veloci- 
ties for a subset of Sample 2; only companions of 
galaxies with —19.5 < Mb < —20.5 are shown. 
Figure 7 shows a comparable plot for Sample 
3, but only for the brighter primaries, since the 
statistics are poor for the fainter primaries. The 
raw r.m.s. velocity is shown by the dashed curve. 
Both figures show that the r.m.s. velocity of the 
raw (uncorrected) data either does not change or 
slightly increases with the distance from the pri- 



Fig. 6. — The line-of-sight velocity differences of 
satellites in SDSS Sample 2 (open circles) for pri- 
maries in the range —19.5 < Mb < —20.5. The 
dashed curve shows the r.m.s. velocities for raw 
data with no removal of interlopers (except for the 
±500 km s _1 cut of the velocities). The r.m.s. ve- 
locity is basically flat; it does not depend on the 
distance to the primary. The full curve shows the- 
oretical prediction for equilibrium NFW halo with 
mass M v i r = 5 x 10 12 M . The decline in velocities 
is not observed in the raw data. 

mary. A few years ago this would have been a 
strong argument in favor of dark matter in the 
outer parts of galaxies, but as discussed in the in- 
troduction, a flat r.m.s. velocity curve extending 
to 300 kpc strongly contradicts any existing dark 
matter model. Taken at face value, the raw r.m.s. 
velocities presented in Figures 6 and 7 contra- 
dict the dark matter scenario, and actually favor 
MOND. 

However, if we remove interlopers using the 
Gauss + constant technique discussed in the pre- 
vious section, as shown by the full lines with error 
bars in Figures 7 and 8, we find that the velocity 
dispersion of satellites declines with distance. We 
note that the results for Sample 1, which extends 
to fainter satellites, are similar to those shown here 
for the other two samples. Using Monte Carlo sim- 
ulations, we reject a hypothesis that the r.m.s. ve- 
locity is constant at the 97% level corresponding 
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Fig. 7. — The same as in Figure 6, but for Sam- 
ple 3 with primaries in the range —20.5 < Mb < 
—21.5. The short dashed curve shows the r.m.s. 
velocities for raw data; the r.m.s. velocity in- 
creases slightly with the distance. The full curve 
with the error bars shows the r.m.s. velocity after 
removal of interlopers in each radial bin; the r.m.s. 
velocity is clearly declining with distance. It is 
consistent with the NFW prediction (full curve) 
for halos with virial mass, M v i r = 6 x 10 12 M©. 
The long dashed curve shows the escape velocity 
from the NFW halo of this mass. 

approximately to the 3cr level for Gaussian statis- 
tics. For the simulations we combine Samples 2 
and 3, which provide independent measurements. 
We assume that deviations of the r.m.s. velocities 
are Gaussian in each of four bins for each sam- 
ple. The average value of the r.m.s. velocity is 
defined by four random points. The probability 
of having \ 2 larger than the observed x 2 and de- 
clining curves in two independent realizations is 
2.9 x 10~ 2 , based on 10,000 realizations. Most of 
the constraints come from Sample 2, where the 
first and the last points deviate from the average 
by 1.5ct and 1.8<r correspondingly. Sample 3 alone 
is consistent with the constant r.m.s. velocity. 
Still, even in this case the velocity declines. When 
we combine both samples, the constant r.m.s. ve- 
locity becomes very improbable. 

The number of estimated interlopers is quite 



Fig. 8. — The same as in Figure 6, but with the 
removal of the interlopers. The full curve with the 
error bars shows the r.m.s. velocity after removal 
of interlopers. As in Figure 7, the r.m.s. is clearly 
declining and consistent with the NFW profile 
(short-dashed curve) with M v ; r = 1.5 x 10 12 M Q . 
The long dashed curve shows the escape velocity 
from the NFW halo of this mass. All satellites 
above the escape velocity curve are interlopers. 

modest. For Sample 1, which has 132 satellites for 
primaries in the range —19.5 <Mb< —20.5, we es- 
timate that 23 are interlopers: 17% of the sample. 
The numbers are similar for the other samples: 
19% in Sample 2 and 20% in sample 3. As ex- 
pected from geometrical arguments, the number of 
interlopers increases with the radius. For example, 
in Sample 3 the number of estimated interlopers 
is 16.7 for radii in the range 100-200 kpc, and 54 
for 200-360 kpc, which is roughly consistent with 
a constant number of interlopers per unit volume. 

The observed decline in velocity dispersion with 
radius closely matches that predicted by an NFW 
dark matter density profile. Figures 7 and 8 give 
the predicted NFW velocities for halos of mass 
6 x 10 12 Mq (for primaries with -20.5 < M B < 
-21.5) and M vir = 1.5 x 10 12 M Q (for primaries 
with -19.5 < M B < -20.5). 

If we remove interlopers using the method of 
McKay et al. (2002), the mean velocity disper- 
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sion for each magnitude bin is smaller than the 
raw value, but we do not get a declining veloc- 
ity curve. This explains the difference between 
our conclusions and those of McKay et al. (2002), 
who do not find a radial dependence on velocity 
dispersion; as discussed above, we feel that our 
procedure of allowing for the radial dependence of 
interloper contamination is better motivated phys- 
ically. 

4.2. Relation of satellite velocity disper- 
sion with host luminosity 

A comparison of Figures 7 and 8 clearly indi- 
cates that the r.m.s. velocity increases with the 
luminosity of the primary. This contradicts the 
earlier results of Zaritsky et al. (1997). Because 
the NFW profile provides an accurate fit to the 
observed velocity dispersion profile, we can use it 
to estimate the virial mass of each group of galax- 
ies. For galaxies with —19.5 < Mb < —20.5 we 
find M v ; r = 1.5 x 10 12 M©. For the brighter mag- 
nitudes, —20.5 < Mb < —21.5, the virial mass is 
M vir = 6x 10 12 M©. These values give us virial 
mass-to-light ratios M/L = 100 and M/L = 150 
for the dimmer and brighter galaxies correspond- 
ingly. Thus, the M/L ratio increases with the lu- 
minosity roughly as M/L L 5 . We believe that 
this change in M/L with luminosity is real and 
does not reflect a change in the stellar populations 
for both magnitude bins. A reason for this is that 
we did not see a dependence of the fraction of E- 
S0 with luminosity for the primaries in Sample 1, 
at least from —19 to —22 blue absolute magni- 
tudes. A similar increase of M/L with luminosity 
(M/L oc £°- 4±0 - 2 ) has been seen from weak lcnsing 
studies by Guzik & Seljak (2002). See, however, 
McKay et al. (2001) who, also from weak lensing, 
found no dependence of the M/L on luminosity. 

Instead of trying to estimate masses of the 
galaxies (which is bound to be model-dependent), 
we can study the more direct dependence of the 
r.m.s. velocity on the luminosity of the primary 
galaxy. This can be viewed as an analog of the 
Tully-Fisher relation. Figure 9 presents the de- 
pendence of the satellite r.m.s. velocity, computed 
within 120 kpc as a function of galaxy luminosity. 
We calculate the r.m.s. velocity relatively close to 
the primary galaxy because within this distance 
the velocity dispersion depends only weakly on 
the distance to the primary for all primaries con- 



sidered here. Hence, we can meaningfully com- 
pare results from galaxies that might have differ- 
ent sizes. The estimated value of the velocity dis- 
persion for an L* galaxy from weak lensing studies 
(Hoeskstra, Yee, & Gladders 2001; McKay et al. 
2001) is also shown; the results of lensing stud- 
ies are very close to those derived from the SDSS 
satellite data. 

The results are consistent with a slope of a cx 
L 0,3 . We can write this dependence using a form 
analogous to the Tully-Fisher relation: 

M B = 0.02-8.71og(2cr). (1) 

This relation is shown in Figure 9. The slope of 
this relation is the same as for spiral galaxies (Ver- 
heijen 2001), which is shown as a long-dashed line 
in the figure. The zero point for the spiral rela- 
tion was shifted down by factor of 1.6 to match the 
SDSS data points; this factor is close to the naively 
expected \/3 correction from rotational velocity to 
ID line-of-sight velocity. The remaining difference 
of 8% is difficult to assign to any particular ef- 
fect; there are several effects that could produce 
the small difference (e.g., projection effects, and 
non-isotropic orbits). 

However, when we compute the r.m.s. velocity 
within a larger radius of 350 kpc, we find a cx L 5 . 
The reason for the steeper slope is a combination 
of the very large radius and declining a(r), which 
will occur if the density profiles of halos at different 
mass are not homologous, e.g., if the concentration 
of the halo decreases with increasing mass. This 
is consistent with theoretical prediction. 

It is rather straightforward to find the cr-virial 
mass relation for the NFW halos. Using the model 
described in Section 3.1 (equilibrium NFW with 
slightly anisotropic velocities), we find the line-of- 
sight r.m.s. velocity averaged for radii 20-100 kpc 
and the line-of-sight velocity for projected distance 
350 kpc. The results are shown in Figure 10. For 
the range of masses considered here, the a-M rela- 
tions are power laws. Just as the observed slope of 
the a-L relation, the slope of the a-M vlv relation 
increases with the radius within which the r.m.s. 
velocity is calculated. For the 350 kpc radius 
the slope 0.5 is the same as the slope of the a-L 
relation for the same radius. This would imply 
a constant mass-to-light ratio. This result, being 
formally correct, is a bit misleading. Measured in 
terms of a more physically motivated quantity, the 
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virial radius, the M/L ratio increases with lumi- 
nosity. 

Comparison with theoretical predictions of the 
luminosities is more complicated and will be de- 
ferred to another paper. We use only the virial 
M/L for two magnitudes Mb = — 20 and —21. In 
a recent paper, Yang et al. (2002) ask which model 
parameters are required for the cosmological mod- 
els to explain the observed luminosity function 
and the Tully-Fisher relation. The main uncer- 
tainty is the M/L ratio. Yang et al. (2002) pre- 
dict an M/L(V) that is needed to account for ob- 
servations. We compare our two M/L ratios with 
the Yang et al. (2002) predictions and find them 
remarkably close. However, this comparison ac- 
tually required numerous corrections, including a 
correction for overdensity of 180 to the virial over- 
density, correction for galactic absorption, correc- 
tion for different bands, and scaling with the Hub- 
ble constant. The corrections were provided by 
Rachel Somervillc. Preliminary comparison with 
semi-analytical models (Somerville, private com- 
munication) indicate a problem: observed galax- 
ies are «0.75 mag too bright as compared with 
theoretical predictions. 

5. Conclusions 

We construct samples of isolated galaxies and 
associated satellites using data from the SDSS. 
We detect about 3000 satellites with absolute blue 
magnitudes going down to Mb = —14. Analy- 
sis of the satellite velocities clearly indicates that 
galaxies are embedded in large halos extending to 
distances up to 350 kpc. Our measured velocity 
dispersion for ~ L* galaxies is in agreement with 
measurements from weak lensing studies. 

Our analysis of galaxies selected from the SDSS 
database clearly indicates that the line-of-sight 
r.m.s. velocity of satellites declines with the dis- 
tance to the primary galaxy. Similar results were 
found for each of our three samples, suggesting 
that it is robust to statistical effects and the ex- 
act definition of an isolated galaxy. This decline 



agrees remarkably well with p cx 



predicted 



by all cosmological models. Both the isothermal 
and MOND profiles contradict the observational 
results. Observations of weak lensing from the 
SDSS combined with the Tully-Fisher and Fun- 
damental Plane relations (Seljak 2002) also finds 




Fig. 9. — The dependence of satellite r.m.s. veloc- 
ity on the absolute B magnitude of the primary 
galaxy. The solid line shows r.m.s. velocity inside 
120 kpc. The long-dashed line shows the Tully- 
Fisher relation of spiral galaxies (Verheijen 2001) 
scaled down by a factor of 1.6 to account for the 
difference in definitions of velocities. The short- 
dashed line shows the r.m.s. velocity for satellites 
in the outer radial shell 250-350 kpc. For these 
large radii the r.m.s. velocity is consistent with 
a oc L 5 (constant M/L). The asterisk shows the 
weak lensing results for galaxies (Hoeskstra, 
Yee, & Gladders 2001) 

that the galactic halo profiles has to be steeper 
than isothermal at large radii. 

Interlopers play an important role. We believe 
that the decline in satellite velocity was not seen 
previously (Zaritsky & White 1994; Zaritsky et al. 
1997; McKay et al. 2002), mainly because of ei- 
ther lack of statistics or insufficient removal of in- 
terlopers. We use numerical simulations to test 
prescriptions for correcting the effects of interlop- 
ers and isolation criteria. Both effects - interlopers 
and isolation criteria - have a tendency to overes- 
timate of the mass of the dark matter halo and 
imply a velocity profile that is too flat. 

We find that the satellite velocity dispersion 
inside a projected 100 kpc radius and the abso- 
lute blue magnitude of the galaxy are related as 
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Fig. 10. — Theoretical predictions for the line-of- 
sight r.m.s. velocity, a, at different projected dis- 
tances as the function of galaxy virial mass. The 
dots show the velocities for equilibrium NFW ha- 
los. The top curves and symbols are for velocities 
averaged for radii 20-100 kpc. The bottom curves 
are for projected distance 350 kpc. Power-law fits 
indicated in the plot with full lines provide very 
accurate approximations. The dashed curves are 
ID tangential r.m.s. velocity at 100 kpc and at 
350 kpc. 

M = 0.02 - 8.71og(2cr) (a cx i 3 ). This depen- 
dence is in agreement with the slope and the zero 
point of the Tully-Fisher relation for spiral disks. 
If we calculate the r.m.s. velocity within a radius 
of 350 kpc, the slope of the velocity-luminosity 
relation is visibly steeper: u cx L 5 . The likely 
explanation for the difference in the slopes is an 
interplay between the shape of the velocity disper- 
sion and the virial radius of galaxies. 

The mass-to-light ratio increases with luminos- 
ity as M/L cx L 5 . For M B = -20 (-21) we find 
that the virial M/L = 100 (150). 
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